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Context, rj Chamaeleontis is a unique young (~ 9 Myr) association with 18 systems concentrated in a radius of x 35 arcmin, i.e. Ipc 
at the cluster distance of 97pc. No other members have been found up to 1.5 degrees from the cluster centre. The cluster mass function 
is consistent with the IMF of other rich young open clusters in the higher mass range but shows a clear deficit of low mass stars and 
brown dwarfs with no objects below 0. IMg. 

Aims. The aim of this paper is to test whether this peculiar mass function could result from dynamical evolution despite the young 
age of the cluster 

Methods. We performed N-body numerical calculations starting with a log-normal IMF and different initial conditions in terms of 
number of systems and cluster radius using the code NB0DY3. We simulated the cluster dynamical evolution over 10 Myr and 
compared the results to the observations. 

Results. We found that it is possible to reproduce r] Cha when starting with a very compact configuration (with A',,,,, = 40 and 
Rq = O.OOSpc) which suggests that the IMF of the association might not be abnormal. The high initial density might also explain the 
deficit of wide binaries that is observed in the cluster 

Key words. Stars : low-mass, brown dwarfs - Stars : mass function - Open clusters and associations : individual : Eta-Chamaeleontis 



> 

oo 

l> 
O 

o 



X 



1. Introduction 

The recently-discovered 77 Chamaeleontis stellar cluster 
(Mamajek et al. 1999) is remarkable in several respects. It is 
a compact (extent ~ Ipc) and sparse group of 18 pre-main 
sequence (PMS) systems, and it is one of the nearest clusters 
to the Sun with an Hipparcos distance d ^ 97pc. The stellar 
members present a high degree of coevality with an age of 
8-10 Myr (Lawson & Feigelson 2001, Luhman & Steeghs 
2004, Lyo et al. 2004a) and were probably born near the Lower 
Centaurus Crux subgroup of the Sco-Cen OB association 
(Mamajek et al. 2000, Jilinski et al. 2005). Despite the youth 
of the cluster there is no evidence for extinction due to remnant 
molecular material, and the proximity of the cluster ensures 
that foreground reddening is negligible (Lyo et al. 2004a). 
Because of these characterictics, the 77 Cha cluster is an ideal 
target to study the properties of "older" T-Tauri stars, such as 
circumstellar disk evolution (Bouwman et al. 2006) and binarity 
(Lyo et al. 2004b, Brandeker et al. 2006), and to search for its 
brown dwarf population. 

Several photometric surveys have been performed during the 
last few years with the aim of obtaining a complete census of 
cluster members beyond the original X-ray-selected population. 
To date, 18 systems have been discovered with a mass range of 
O.lSMo to 3.8M0 (Mamajek et al. 1999, Lawson et al. 2002, Lyo 
et al. 2004b, Song et al. 2004). Lyo et al. (2004b) found that the 
cluster mass function is consistent with the IMF of rich young 
open clusters and field stars in this mass domain. By extrapolat- 
ing the mass distribution of the Trapezium cluster (Muench et 
al. 2002) to lower masses, they predicted 10 - 14 additional low 
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mass stars with O.O8M0 < m < 0.15Mq and about 15 brown 
dwarfs should accompany the documented stellar population. 
However, Luhman (2004) examined the region within a radius 
of 1.5° surrounding t] Cha using 2MASS and DENIS photo- 
metric data and found no new members within the mass range 
0.025 - O.IOMq. Lyo et al. (2006) performed a deeper survey 
of the central core region down to ~ 13My„y, and found no new 
low mass members either. Whether this clear deficit of very low 
mass stars and brown dwarfs in the cluster mass function results 
from an abnormal IMF or is due to dynamical evolution is not 
completely clear yet. Although 77 Cha is young and sparse, there 
is some hint that dynamical evolution has already occurred. In 
particular, the radial distribution of the cluster members indi- 
cates that significant mass segregation is present with more than 
50% of the mass residing in the inner 6 arcmin or 0.17 pc (Lyo 
et al. 2004b). This is the sparsest stellar cluster for which mass 
segregation has been seen. Moreover, the maximum stellar mass 
- cluster mass relation proposed by Weidner & Kroupa (2006) 
yields a mass of 30 - 40Mq for the embedded cluster at the ori- 
gin of 77 Cha. The actual cluster mass is I8/W0, which suggests 
that a significant fraction of the primordial stellar members has 
already been lost, probably as a consequence of gas expulsion. 

Binarity studies indicate that the 18 stellar systems of rj Cha 
include five confirmed (RECX 1, RECX 7, RECX 9, RECX 12, 
and RS Cha which may be a triple) and three probable (ECHA 
J0836.2-7908, ECHA J0838.9-7916, and 77 Cha) binaries; see 
Table 2 of Lyo et al. (2004b). This yields a multiple fraction 
of 28 - 44% which is 2 - 4 times higher than that of solar- 
type field stars, although it is of the same order as observed in 
other young, nearby associations such as Taurus or TW Hydrae. 
Speckle observations (Kohler & Petr-Gotzens 2002) and adap- 
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tive optics imaging (Brandeker et al. 2006) resolved only two 
binaries, RECX 1 and RECX 9. In both cases the projected sep- 
aration is ~ 0.2" which corresponds to about 20 AU at a distance 
of 97 pc, and the other binaries are not spatially resolved. RECX 
7 is a dual-lined spectroscopic binary (Lyo et al. 2003) and RS 
Cha AB is an eclipsing system (Andersen 1991, Mamajek et al. 
2000). No companion with projected separation larger than 20 
AU has been found around the 18 cluster members. Brandeker 
et al. (2006) computed a companion probability based on their 
contrast sensitivity limit and found an upper limit of about 18% 
for wide (> 30 AU) binaries in 77 Cha. This contrasts to the ~58% 
in the TW Hydrae association that is of similar age. Bouwman et 
al. (2006) presents the most recent summary of binary properties 
for the cluster. 

These properties - the deficit of very low mass objects, mass 
segregation, and the lack of wide binaries - are rather unusual 
considering the youth and sparsity of the rj Cha cluster. Whether 
they were imprinted during the formation phase as a result of 
different initial conditions, or whether they are due to dynamical 
evolution is an open issue. 

The aim of this paper is to test this second possibility. We 
try to reproduce the observed 77 Cha properties using N-body nu- 
merical simulations of the cluster dynamical evolution, starting 
with standard initial conditions. In the next section, we explain 
the model we used for the calculations. The results and analysis 
are presented in Section 3, and the conclusion is given in Section 
4. 

2. Numerical simulations 

2.1. IMF 

There is growing evidence that the IMF is independent of ini- 
tial conditions, from the stellar domain down to the substellar 
regime (e.g. Kroupa 2001, Chabrier 2003, Moraux et al. 2007). 
It would therefore be rather surprising to find such a different 
initial mass spectrum for t] Cha with no members with masses 
below 0.1 SMq. Following the log-normal IMF from Chabrier 
(2003) 

(logm - log nio)^^ 
2cr2 



^(m) oc exp - 



where otq = 0.25Mq and cr - 0.52, we estimate the number 
of systems necessary to reproduce the four objects with a mass 
between 1 and 4Mq observed in the clustefl We found that about 
50 systems with masses between 0.01 and 4Mq are mandatory. 
We started a first set of numerical calculations with Njmt - 50 
and different values of initial cluster radius Rq. Then we fixed 
Rq around the values that give the best results in reproducing the 
observations and we performed a second set of simulations with 
Ni„i, = 30, 40, 60, 70 to test the effect of the initial number of 
systems on the results. 

2.2. Initial cluster radius 

Recently, Adams et al. (2006) found a correlation between the 
number of stars in a star forming region and the core radius Rc 
using a compilation of data on embedded clusters from Lada & 
Lada (2003) and Carpenter (2000). Clusters identified by Spitzer 
using infrared-excess sources also follow this relationship (Allen 



' The stars/systems with m > IMq are r] Cha, the RS Cha binary, HD 
75505, and at least one of the three K-stars, RECX 1, 7 and 11. In the 
absence of dynamical masses for any of the K-type stars, their adopted 
masses are somewhat model-dependent. 



et al. 2006). This corresponds to an average surface density of 
cluster members NIR^ that is nearly constant and gives R^ - 
0.3-1.0 pcforA^ = 50 stars. 

In the calculations, the cluster size is defined by the virial 
radius 

Ro = - 

4E 

where Mc is the total mass of the cluster, and E is the total en- 
ergy of the system. To maximize the density, and thus dynami- 
cal interactions, between cluster members, we first adopted the 
value Rq = 0.3 pc as an initial condition. As 77 Cha appears to 
be already evolved dynamically despite its youth - all the low 
mass objects with m < O.IMq have already been lost - we also 
started our calculations with more compact configurations down 
to Ro = 0.005 pc. We chose Rq = 0.3, 0.1, 0.05, 0.03, 0.01 and 
0.005 pc. 

2.3. Initial density distribution 

Significant mass segregation is already present in 77 Cha with 
more than 50% of the stellar mass residing in the inner 0.17 
pc (Lyo et al. 2004b). Similarly, studies of young rich stellar 
clusters (e.g. in the ONC; Hillenbrand & Hartmann 1998) have 
shown that the most massive stars are concentrated in the central 
core of the cluster Bonnell & Davies (1998) have suggested that 
mass segregation may reflect the initial condition of the cluster, 
and is not a consequence of dynamical interactions. To test this 
hypothesis, we started some of our simulations with primordial 
mass segregation, and some without. 

In all cases, we initially select stellar positions so as to follow 
a Plummer model 



p(r) = 



1 



47tRI, [1 +(r/;?p/)2]5/2 



pi 



with Rpi = (37r/16)/?() (Aarseth et al. 1974). To add primordial 
mass segregation, we simply rejected stars with m > 0.5Mq that 
were assigned to r > O.IRq and repeated the random drawing 
for this cluster member. The velocities are scaled so as to en- 
sure global virial equilibrium, although we note that in the mass 
segregated case the initial conditions do not correspond to a so- 
lution of the steady state collisionless Boltzmann equation. We 
note that rapid collisional evolution of this small N system ren- 
ders the question of whether the initial conditions are in precise 
equilibrium of little importance. 

2.4. Initial speeds 

Initially the cluster is assumed to be in virial equilibrium, which 
is characterized by the ratio of kinetic energy to potential energy 
Q = 0.5. The velocity distribution is everywhere isotropic, and 
the generating function corresponding to a Plummer model is 
given by Aarseth et al. (1974). 

2.5. Gas 

Although the details of gas removal are not fully understood, ob- 
servations indicate that clusters older than 5 Myr are rarely as- 
sociated with molecular gas. This is the case of 77 Cha. Mamajek 
et al. (2000) suggested that the superbubble formed by Sco-Cen 
OB winds and supernova remnants terminated star formation in 
the cluster and dispersed its natal molecular gas. This could have 
occured as recently as only a few Myr ago. 
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The gas is usually represented as an external potential in 
cluster evolution models, and is removed a few Myr after the 
star formation process has begun on a timescale of the order of a 
crossing time or less (e.g. Geyer & Burkert 2001; Kroupa et al. 
2001; Kroupa & Bouvier 2003; Bastian & Goodwin 2006). At 
the end of the gas expulsion, the cluster is out of virial equilib- 
rium and the subsequent evolution is a relaxation of the cluster 
as it attempts to return to virial equilibrium. This phase is found 
to last for about 20 initial crossing times tcr (Geyer & Burkert 
2001; Bastian & Goodwin 2006) but we expect it to be short 
(~ 0.1 - 1 Myr) in the case of tj Cha. In order to reproduce the 
observations, the initial configuration should be very compact, 
which suggests that f^r will be very small initially. 

In this study, we consider that the gas has already been re- 
moved and that the cluster has returned to virial equilibrium 
when we start our calculations. This means we start the models 
when rj Cha is already a few Myr old. The dynamical evolution 
of the cluster is then purely due to two-body interactions and that 
is simulated by our N-body numerical calculations. Therefore, 
the models do not have to run for 10 Myr (the maximum esti- 
mated age of T] Cha) but can be as short as a few Myr. 

2.6. Binaries 

No binaries are present in our calculations, even though there are 
eight known or probable binary systems in the rj Cha cluster In 
order to be able to compare the observations to the simulations, 
the masses of the binary components are combined to make a 
single star The simulations thus contain the correct mass in stars 
but omit interactions due to single-binary and binary-binary scat- 
tering. In a forthcoming paper, we shall include an input binary 
population and undertake a detailed comparison with observa- 
tions of the binary systems in tj Cha. 

2.7. Calculations 

For each set of initial conditions summarized in Table [1] we 
made 100 realisations running an adapted version of the code 
NBODY3 written by Aarseth (1999). It is a direct, high precision 
N-body integrator that uses the sophisticated Kustaanheimo- 
Stiefel (1965) two-body regularization to treat close encounters. 
There is no softening parameter. In total, 1900 simulations have 
been performed. 

The stellar system containing initially A^,„„ sources in the 
mass range 0.01 - 4-Mq is placed in the Galactic potential which 
defines a tidal radius 



\4A(A-B)j 

where A and B are the Oort's constants, and Mc is the clus- 
ter mass. When the distance of an object to the cluster centre 
is larger than 2 x r,, it is considered to have escaped from the 
cluster and it is removed from the calculation. Each calculation 
proceeds until the number of bound systems falls below 14 or 
until the age becomes larger than 10 Myr. The masses, positions 
and velocities of the systems which remain bound to the cluster 
are output every ~ 0.1 Myr. There are about 100 outputs per sim- 
ulation on average, i.e. about 190000 output files were produced 
in total. 

3. Results and analysis 

The aim of the N-body simulations is to find initial conditions 
which are able to reproduce the actual state of the i] Cha clus- 



ter, when starting with a number of stellar systems A^,-,,,, imposed 
by a log-normal IMF. Each output at each time step is analysed 
and compared to the observations and several criteria have to be 
fulfilled before it is considered to be a good match. 

3.1. Number of systems 

First, we look at the number of bound systems A^i^^. located 
within 1 pc of the cluster centre in each output. If the number 
is larger than 22 or smaller than 14, we do not consider the re- 
sult as acceptable. 

We used this simple criterion alone for the first set of sim- 
ulations starting with A^,„„ = 50 to estimate the initial cluster 
radius Rq and the initial density required to lose the majority of 
the cluster members through dynamical interactions in less than 
10 Myr. With Rq = 0.3 pc, the cluster is too diffuse and the two- 
body encounters are not effective enough to eject members. Even 
after 10 Myr (the maximum duration of the simulations), more 
than half the systems (> 25/50) remained inside a radius of 1 pc 
for all of the realisations (see Fig.[T]l. 

The situation improves when we start the calculations with 
a more compact configuration. As expected, more objects are 
found outside Ipc and the cluster population spreads out more 
quickly. With A^,„„ - 50 and Ro = 0.1 pc initially, 21 realisa- 
tions out of 100 end up with Nipc < 22 but none of them reach 
Nij„ = 14 in less than 10 Myr. For R^ = 0.01 pc and 0.005 
pc however, most of the calculations are stopped because they 
get to the final state with Nipc = 14 before 10 Myr. We can 
see from Figure |2] that the more compact the configuration, the 
more quickly Nipc - 14 is attained. Bearing in mind that i] Cha 
has an age of ~ 10 Myr and that the cluster lost its gas com- 
ponent probably after a few Myr (corresponding to the time at 
which our N-body calculations start), the most favorable age at 
the end of our simulations to reproduce t] Cha should roughly be 
between 4 and 8 Myr The results from Figure|2]suggest that the 
best initial cluster radius Ro is between 0.01 and 0.005 pc. We 
thus eliminated all the other values for Ro as initial conditions 
and we retained Ro - 0.005 and 0.01 pc to run the second set 
of simulations with A^,„,-, - 30, 40, 50, 60 or 70 in order to find 
the best value of A^,„,-, capable of matching the observed cluster 
characteristics. 

By applying the same test on the results of these simulations, 
we found that most of them reach 14 < A^i^c < 22 before 10 Myr 
apart from Ni„i, = 70 and Ro = 0.01 pc where only 66 calcula- 
tions fulfill this criterion (see Table |2]i. This is due to the fact 
that there are more objects in the cluster and that it takes longer 
to remove most of them from the centre. This indicates that ini- 
tial conditions with A^,„„ > 70 are not favorable to reproduce i] 
Cha. Then we looked at the age distribution of the calculations 
which ended with 14 < Nipc < 22 (Figure |3j. We found that 
most clusters starting with A^,„,-, - 30 lost their members very 
quickly. There is a broad peak around ~ 2 - 3 Myr in the dis- 
tribution, which is probably too young to reproduce t] Cha. This 
suggests that simulations starting with A^,„„ < 30 are less likely 
to match the observations. However, all of the other initial con- 
ditions with Ni„i, - 40, 50 or 60 and Ro - 0.01 or 0.005 pc seem 
to be a good match to the observed properties according to the 
first criterion based solely on the number of objects left inside a 
radius of 1 pc. 
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Table 1. Summary of initial conditions for each simulation 
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Fig. 1. Number of systems inside a radius of 1 pc for Nmn = 50 and Ro - 0.3, 0. 1, 0.05, 0.03, 0.01 and 0.005 pc at the end of the 
calculation. This corresponds to an age of 10 Myr or less if the number of systems in the simulation falls below 14. 
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Table 2. Number of realisations out of 100 yielding Nipc < 22 
within 10 Myr for N^, = 30, 40, 50, 60, 70 and = 0.01 and 
0.005 pc. 



3.2. Number of Stars withm > IMq 

When starting with a compact configuration (Rq < 0.01 pc), the 
two-body interactions are strong enough to eject a large number 
of objects in a few Myr so that less than 22 stars remain inside a 
radius of 1 pc from the cluster centre. However, we do not want 
this to be the case for the most massive stars as four objects with 
a mass between 1 and 4Mq are observed in the cluster. Therefore 
the second criterion we used is to check that at least 3 stars more 
massive than I.5M0 are included in N]pc, i.e. N„,>i,5Ma ^ 3. 

This second test allows us to eliminate more than half of the 
calculations that fulfill the first criterion, but does not really help 
us to define the best initial conditions as the trend is the same for 
all of them. In particular, the proportion of realisations which 



pass this test for a given A^,„„ is similar for Rq - 0.005 or 0.01 
pc (see TableO. 
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Table 3. Number of realisations out of 100 ending up with 
Nipc < 22 and A^,„>i.5Mo > 3 for Nm, = 30, 40, 50, 60, 70 and 
Ro = 0.01 and 0.005 pc. 



3.3. Number of remaining brown dwarfs 

The third criterion is the comparison between the lower mass 
function resulting from the simulations and the observations. 
No brown dwarf or very low mass (VLM) star with m < 
O.IMq has been found within 1.5 degrees of the cluster centre 
(Luhman 2004) and this has to be reproduced by the simula- 
tions. Therefore, we require that Nvlm, the number of objects 
with m < O.IM0 inside 1.5 degrees, is or 1 to consider a result 
as acceptable. 
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Fig. 2. "Age" at the end of the simulations started with A^,„„ - 
50 and Ro = 0.03 (dotted hne), 0.01 (solid line) and 0.005 pc 
(dashed line). Here "age" means time at which A^i^c dropped to 
< 14 (if this is < 10 Myr) or 10 Myr (if N\pc < 22 at an age of 
10 Myr). The remaining simulations (for which A^i^c > 22 at 10 
Myr) are not shown here. 
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Fig. 3. Age at the end of the simulations started with A^,„„ - 30, 
40, 60 and 70 which end up with 14 < N\pc < 22. In each panel 
the solid histogram corresponds to Rq - 0.01 pc and the dotted 
one to 0.005 pc. 

This is a very stringent constraint. We saw in Section 13.11 
that it is not too hard to eject a large number of cluster members 
from the centre when starting with a very compact configuration. 
Almost all the realisations with Ro - 0.01 or 0.005 pc end up 
with 14 < N\pc < 22. It is already more difficult to keep the most 



massive stars in the core - only half of the reaUsations passed this 
test - but it is even harder to get all the very low mass objects far 
away from the cluster centre in less than 10 Myr. Table |4]gives 
the number of realisations which fulfill criteria 1 and 3 alone 
(14 < Nypc < 22 and Nvlm < 1) and all of criteria 1,2 and 3 
(14 < N\pc < 22, Nm>i.5Ms ^ 3 and Nvlm ^ 1)- We can see 
that the number drops down, around 10 or even less for some 
initial conditions, even if we relax the second test on N„,>i,5Mo- 
This suggests in particular that Ro - 0.01 pc is not favorable to 
reproduce 77 Cha, apart perhaps for A^,„/, - 40. 
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Table 4. Number of realisations out of 100 ending with 14 < 
Nipc < 22 and Nvlm ^ 1 (first column) and with 14 < Nipc < 
22, Nm>\.5Mo ^ 3 andA^vLM ^ 1 (second column) for each initial 
conditions Nm, = 30, 40, 50, 60, 70 and Ro = 0.01 pc, 0.005 pc 
and 0.01 pc with primordial mass segregation. 



The reason why it is so difficult to fulfill the three criteria 
at the same time is that the loss rate becomes roughly the same 
for all stars, independently of their mass, when the initial con- 
figuration is very compact. In order to keep the massive stars in 
the centre, the dynamical interactions should not be too strong to 
kick them out but then, even if some brown dwarfs escaped the 
cluster core, several still remain. In less than 10 Myr, the mass 
segregation favoring the evaporation of the lower mass objects 
does not really have time to occur. Simulations starting with a 
large number of systems A^,„„ like 60 or 70 have a larger number 
of brown dwarfs initially and it is very difficult to move them 
away in less than 10 Myr. Starting with Rq = 0.05 pc instead of 
0.01 pc helps to some extent as the configuration does not get 
too compact and the massive star are not removed from the core 
(which is the case for Ni„i, - 70 and Rq = 0.005 pc). 

In order to reproduce the mass distribution of rj Cha, one 
needs to lose more brown dwarfs and VLM stars than higher 
mass stars. We can therefore wonder if primordial mass segre- 
gation would help as it could yield such a preferential escape 
rate. However, this does not seem to be the case: the number of 
good realisations with R^ - 0.01 pc does not really get larger 
with primordial mass segregation than without (see TableHJ. In 
some cases (A^,-,,;, = 30 and 40) it even gets smaller This suggests 
that primordial mass segregation is not required to reproduce 77 
Cha. In order to check this result, we performed a last test on 
our simulations concerning the radial distribution of the cluster 
members. 

3.4. Radial distribution 

We searched the ROSAT All-Sky Survey (RASS) for cluster 
members with projected separations of up to 5 degrees. For 10 
Myr-old stars T Tauri stars at ~ 100 pc, the RASS is sensitive 
to mid-K to early-M stars, i.e. to near-solar mass stars. To be 
candidate members, RASS-selected objects should show high 
levels of lithium indicative of young solar-type stars, photom- 
etry in agreement with the placement of known cluster members 
in colour-magnitude diagrams, and three-dimensional space mo- 
tions consistent with ejection from the cluster in the past 2-10 
Myr, at a velocity less than ~ 5 km/s. Only rarely is this set 
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of quantities complete and available with precision, e.g. optical 
photometry derived from plate scans is poor and radial veloc- 
ities are usually absent. Proper motions are available from the 
UCAC2 catalogue, while a high-resolution spectroscopic survey 
of RASS sources across the Chamaeleon region was performed 
by Covino et al. (1997) which included a few known 77 Cha 
stars. Within proper motion vectors < 10 mas/yr of the cluster's 
proper motion of = (-29.9,27.5) mas/yr (Mamajek et 

al. 1999), we found only a few additional RASS stars that might 
meet these criteria for membership. Stars ejected from the 77 Cha 
with a velocity higher than 10 mas/yr, corresponding to 5 km/s 
at the distance of the cluster, will therefore not be sampled but 
our simulations show that they are rare. High velocity objects are 
ejected in the early stage of the cluster evolution and are already 
at a much larger distance than 5 degrees from the cluster center. 
While noting that these surveys are imprecise and/or incomplete, 
this outcome indicates there is not a large number of undiscov- 
ered solar-mass members within 5 degrees of the cluster core. 

All of the cluster members are instead concentrated within 1 
pc from the cluster centre and only a few additional solar type 
stars, if any, can be located within 5 degrees. Therefore we re- 
quire that the number of stars with m > O.SM© located between 
1 and 8.75 pc (= 5 degrees for a distance of 100 pc) from the 
cluster centre is less than 2. This constitutes our last criterion 
to check the ability of the simulations to reproduce the observa- 
tions. 

When performing this test on the realisations that already 
fulfill the three previous criteria, the number of remaining good 
results becomes very small for all the initial conditions (see 
Table|5]l. Again we find that starting with primordial mass segre- 
gation can help in some cases but it is not necessary. A few per- 
cent of the simulations starting with Rq - 0.005 pc can always 
reproduce 77 Cha independently of A^,„„ which suggests that this 
corresponds to the best initial condition for the cluster radius. 
Similarly, A^,„/, - 40 is the most favorable initial number of sys- 
tems, as we obtain good results for both R() - 0.005 and 0.01 pc. 
An example of a good realisation resulting from a calculation 
started with A^,„» = 40 and Rq - 0.005 is shown in Figure|4] 
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Table 5. Number of realisations out of 100 fulfilling the full set 
of selection criteria for each initial conditions A^,-,,,, - 30, 40, 
50, 60, 70 and Rq - 0.01 pc without and with primordial mass 
segregation and Rq = 0.005 pc. 



It is important to note, that even if only a few percent of the 
realisations fulfill all the above criteria, they remain valid across 
several time steps, i.e. the time for which we have a good match 
with the observations last for more than 1 Myr, typically at an 
age between 4 and 10 Myr. For example, in the realisation from 
which Fig. |4]is drawn, the simulated cluster remains in agree- 
ment with the observations between 3.5 and 10 Myr The fact 
that such a dynamical state is not transient indicates that the cal- 
culations can mimic an observational state and that our simula- 
tions can effectively reproduce t] Cha. The difference between 
the age given by the simulation and the cluster age would then 
correspond to the time at which the gas was lost. 



4. Conclusion and discussion 

4.1. Summary of results 

To summarize, the criteria we used to consider that a result is 
consistent with the observations are the following: 

- the number of objects within 1 pc of the cluster centre is 
14 < Nipc < 22 at some point at f < 10 Myr, 

- it includes at least 3 stars more massive than 1 .5Mq, 

- there is no more than one brown dwarf or VLM star with 
m < 0.1 M0 within 1.5 degrees, 

- and less than 2 stars with m > Q.5Mq are located between 1 
and 8.75 pc (corresponding to 5 degrees at a distance of 100 
pc). 

We tested initial conditions with Ro - 0.3 pc down to 0.005 
pc and Ninit = 30, 40, 50, 60 and 70 distributed according to a 
log-normal IMF. None of the simulations starting with Ro > 0.01 
pc give good results, except if Rq = 0.01 pc and A^,„„ - 40. For 
the calculations with Rq = 0.005 pc however, there are always a 
few realisations that match the observations whatever the value 

of Mini,. 

We thus found that it is effectively possible to reproduce t] 
Cha, which suggests in particular that the deficit of brown dwarfs 
and very low mass stars observed in the present day cluster mass 
function may not be due to a peculiar IMF but to dynamical evo- 
lution. 

4.2. Compactness of the core 

The most favorable initial condition in reproducing 77 Cha is 
Ninir = 40 and Rq = 0.005 pc, which corresponds to an ini- 
tial central density of 3.7 x 10'^ stars/pc^. In the example shown 
FiglH the total cluster mass at the beginning of the calculation 
is 21.4Mo. This corresponds to an initial crossing time of 0.003 
Myr which confirms a posteriori our assumption that the relax- 
ation phase after the gas expulsion is very short (< 0.1 Myr). 
According to the maximum stellar mass - cluster mass relation 
discussed by Weidner & Kroupa (2006), the mass of the embed- 
ded cluster at the origin of 77 Cha was of the order of 30 - 4OM0. 
This suggests that 30 - 40% of the stars had been lost when the 
simulation starts and that the initial density may have been even 
higher than the value we obtain. Further calculations including 
the treatment of gas disruption (as well as an input binary pop- 
ulation) will be done in a forthcoming paper to investigate this 
effect. 

Nevertheless, a density of ~ 10^ stars/pc^ is extraordinar- 
ily high and potentially makes 77 Cha a unique region where 
star formation would have occurred in a very compact config- 
uration. It is 5 orders of magnitude larger than the initial den- 
sity taken by Kroupa & Bouvier (2003) to simulate the aggre- 
gates in Taurus and 2-3 orders of magnitude larger than that 
used for the ONC and Pleiades models (Kroupa et al. 2001). In 
the simulations of massive star cluster formation performed by 
Bonnell et al. (2003), the stellar density attains a maximum of 
10^ to 10^ stars/pc^ but this occurs only locally. At the cluster's 
scale the density is typically 10^ stars/pc^, 3 orders of magni- 
tude smaller than what we have. Observationally, such a high 
density has never been found in any embedded cluster but this 
may be due to a resolution effect. An initial radius Rq = 0.005 
pc corresponds to a half mass radius of 0.0035 pc or 7.4" on the 
sky assuming a cluster distance of 97 pc. At the distances to the 
youngest large-scale regions of star formation such as the Carina 
nebula at d - 2.3 kpc, such a cluster would not be spatially re- 
solved (or only marginally) by radio and IR observations and 
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Fig. 4. Example of the result of a realisation considered as acceptable. The initial conditions are A^,„,-, - 40 and Rq - 0.005 pc. Top 
left panel: cluster mass function from the simulation (solid histogram) and the observations (dashed histogram). Top right panel: 
Cluster member spatial distribution. The dashed circles centreed on the cluster correspond to a radius of 1 pc and 2.6 pc (=1.5 
degrees at 100 pc). The stars with m > O.IM© are shown as circles and the symbols are larger for larger mass. The very low mass 
stars and brown dwarfs are shown as crosses. Bottom left panel: Cumulative distribution of the number of cluster members for the 
simulation (solid line) and the observations (dashed line). Bottom right panel: Mass cumulative distribution for the simulation (solid 
line) and the observations (dashed line). 



would instead be seen as a bright single source. High angular 
resolution observations (e.g. with ALMA) are needed to know if 
such clusters exist. 

4.3. Binarity 

For the large majority of systems, the smallest nearest-neighbour 
distance is obtained at f = and is around 5 x 10"^ - lO^^* pc, 
which corresponds to 100 - 200 AU. This is consistent with the 
typical ejection velocities of a few km/s (Fig. |5]) and suggests 
that binaries cannot have separations larger than a factor of a 



few times less than this distance. Kroupa & Burkert (2001) per- 
formed N-body calculations of extremely compact clusters, with 
a central density similar to what we find for rj Cha. They start 
with a binary proportion fh = l and a uniform logarithmic pe- 
riod distribution in the range 4.5 < log P < 5.5, with P in days. 
This corresponds to a separation interval ~ 20 - 100 AU for 
an equal mass binary of a total mass IMq. For an initial den- 
sity logpc - 8.9 stars/pc^^, the binary proportion in this period 
range falls down to - 0.2 very quickly, in a timescale of the 
order of a crossing time. This is consistent with the results from 
Brandeker et al. (2006) who found an upper limit of 1 8 % of wide 
binaries at projected separations > 30 AU, and also with the ap- 
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parent absetice of circumbinary disks in Spitzer IRS spectra of 
Cha (Bouwman et al. 2006). 

4.4. Where are the escapers ? 

The objects that reach twice the tidal radius, i.e. r = 6 to 8 pc de- 
pending on the simulations, escape the cluster and are removed 
from the calculations. But it is still possible to follow their evolu- 
tion afterwards because we know their position, velocity and the 
time at which they escape. We simply assume that their move- 
ment is linear The velocity of the escapers is decreasing with 
time as the density gets smaller and the interactions are softer 
We find that it goes roughly with 1 /time and that it is indepen- 
dent of the object mass, i.e. the ratio of VLM objects to normal 
stars is roughly constant with radius and is consistent with the 
input IMF. All the systems that are ejected in the early phase 
(the first Myr) have a velocity of ~ 6 km/s or more which means 
that they are now more than ~60 pc (or 30 degrees) away, given 
the age and distance of 77 Cha. Most of the objects are ejected be- 
tween 1 and 4 Myr with a velocity in the range 1-5 km/s. If we 
assume that the simulations run for 7 Myr to reproduce the clus- 
ter as in the example shown in Fig.lH they will be at a distance 
of at least 9 pc. The systems that are ejected later have a smaller 
velocity and remain closer to the cluster centre as there is less 
time to spread them away. Statistically, for all the simulations 
starting with A^,„„ - 40 and Rq = 0.005 pc, we find a larger num- 
ber of systems around 6 - 10 pc in projected distance (or 3.5 - 6 
degrees for a cluster distance of 97 pc) whatever the time is (see 
Figure |5]l. The surface density is so small (~ 0.02 objects/pc^) 
that it would be quite difficult to observe this population pho- 
tometrically. However proper motion surveys might be able to 
select for ejected cluster members as their projected velocities 
should be distributed around the cluster motion of (-29.9, 27.5) 
mas/yr, which is quite high. Recently Cruz et al. (2007) discov- 
ered a population of young (5 - 50 Myr) brown dwarfs in the 
field. A few of them could come from 77 Cha and would have 
been ejected early with a high velocity (> 6 km/s). 
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Fig. 5. Statistics over all the 100 simulations starting with A^,„„ = 40 and Rq = 0.005 pc. Top left panel: Average number of 
escaped systems per time bin (0.5 Myr) with m > O.IMq (solid line) and m < O.IM© (dashed line) as a function of time. Top right 
panel: Ejection velocity of all the escapers versus time (for all the 100 realisations). Crosses are VLM stars and brown dwarfs with 
m < O.IMq. Bottom left panel: Radial distribution of the escapers at t-7 Myr for m > O.IMq (solid line) and m < O.IM© (dashed 
line). The number of objects is averaged over the 100 realisations and is given per radius bin of 1 pc. Bottom right panel: Proper 
motion diagram of all the escapers at t=7 Myr. 



